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ABSTRACT 

We have performed 2D hydromagnetic simulations with an adaptive mesh refinement 
code to examine the response of a pre-existing initially spherical dense core to a non- 
linear fast-mode wave. One key parameter is the ratio of the wavelength to the initial 
core radius. If that ratio is large and the wave amplitude is suffic i ent, s ignificant 
compression of the core occurs, as envisaged by iMvers fc Lazaria"nl (|1998[ ) in their 
"turbulent cooling flow" picture. For smaller values of that ratio, an initial value of 
the ratio of the thermal pressure to magnetic pressure of 0.2, and sufficiently large wave 
amplitude, the scattering induces the production of dense substructure in the core. 
This substructure may be related to that detected in the dense core associated with 
the cyanopolyyne peak in TMC-1. Our simulations also show that short- wavelength 
waves, contrary to large- wavelength waves, do not confine dense cores. 

Key words: MHD - stars: formation - ISM: clouds - ISM: individual: TMC-1 - 
ISM: molecules 



1 INTRODUCTION 



iBallesteros-Paredes fc MacLowl (l2002MPadoan fc Nordlundl 



d2002h, iGanimie et all (|2QQ3l), iLi et alJ (120041)," 
INakamura fc Lil (l2005h . and I Vazquez- Semaden i et al.l 



20051 ) are amongst the authors who have performed 3D 
simulations to investigate the formation of clumps in media 
in which the value of /3, the ratio of the thermal pressure 
to magnetic pressure, is initially small everywhere but in 
which velocity perturbations exist. These calculations have 
included many ingredients. 

In th ree previous pa p ers, we have carried 
out ID (iFalle fc Hartauistl l200l hereafter FH02; 
iLim Falle fc Hartauistl I 20051. herea fter LFH05) and 
2D (|Van Loo. Falle fc Hartauistl 120061 . hereafter VFH06) 



numerical studies of MHD waves in low-/3 plasma to elu- 
cidate the mechanism by which clumps form. We consider 
the results of pa rticular relevance to the formation of 
dense cores (e.g. ICaselli et al. I l2002h within translucent 
clumps like thos e discovered in CO maps of the Rosette 
Molecular Cloud (jWilliams. Blitz fc St arkll 19951 ). The FH02 
and VFH06 studies concerned small-amplitude waves and 
involved ID and 2D calculations, respectively. A significant 
difference in the results of the two studies is that dense 
core-like structures persisted longer in the 2D simulations 
than in the ID simulations. In 2D, such a structure 
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survived as long as the amplitude that might be loosely 
associated with the fast-mode component of the large-scale 
perturbation remained a significant fraction of its initial 
value. 

The interaction of fast-mode waves with large-scale 
density inhomogeneities produced sub-structure within the 
large-scale clumps. The large-scale structures themselves 
arise due to the non-linear steepening of the fast-mode 
waves. The present paper is concerned with a related but 
somewhat different problem: the response of a pre-existing 
dense core-like structure to a background wave. The effect 
of large-amplitude hydromagn etic waves has alr eady been 
studied analytically by iMcKee fc Zweibell ([l992). Using a 
time-averaged form of the virial theorem, they showed that 
a turbulent medium can pressure- confine a dense clump. A 
limitation of their analysis is that it does not describe the 
effect of large-wavelength perturbations. With our numeri- 
cal model, we can study the effect of both large and small 
wavelength perturbations. 

Numerical si mulations of this problem are li mited, as 
most studies (e.g. iMac Low. McKee fc Kleinl lT994) only de- 
scribe the effect of a shock hitting an inhomogeneity. The 
flow initially compresses the dense core, but eventually de- 
stroys and fragments the core. A similar behaviour can also 
be expected for non-linear waves. Our simulations can be 
considered a natural extension and generalisation of work 
on shocks interacting with inhomogeneities. 

The response of a dense core to a background wave may 
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well be relevant to the detection of substructure in Core D, 
w hich coincides with the cyanopolyyne peak, in the TMC-1 
bv lPeng et all (jl998l ). TMC-1 is a ridge of dense cores, each 
of which has a size of 0.1 - 0.2 pc and a m ass of one to several 
solar masses (e.g. iHirahara et al.l[l992h . The substructures 
studied by IPeng et al.1 (|l998l ) have sizes of 0.03 - 0.06 pc 
and masses of 0.01 - 0.15 Mq. Most are not bound by self- 
gravity and are too small to produce brown dwarfs. Since 
many of these substructures have masses below the Jeans 
mass, they must have formed by some mechanism other than 
gravitational instability. 

A possible mechanism for forming the substructure may 
be the same as or related to the mechanism giving rise 
to dense core formation in the picture that we have pur- 
sued in FH02, LFH05 and VFH06 and which is, no doubt, 
contributing to dense core formation in the other numeri- 
cal simulations referenced above. However, that mechanism 
ceases to give large density contrast structures if (3 ^ 1 un- 
less the amplitude of the wave is hyper- Alfvenic (FH02, 
LFH02). So me dense cores may have v a lues of (3 as low 
as 1 (e.g. IWard- Thompson et al.1 l200d . I Ward- Thompson! 
20021) but some may have values of 3 greater than unity (e.g. 
Kirk. Ward- Thompson fe Crutch er 2006). The lines emitted 
by species such as CS, NH3 and CO in dense cores show sub- 
stantial non-thermal components to their broadening, which 
are, however, usually not highly sup ersonic with respect t o 
the thermal sound speed in H2 (e.g. iFuller fc Mversl [ |l992h . 
For the substructure in Core D of TMC-1, the non-thermal 
line broadening component is important but in most sub- 
structures markedly subsonic in the above sense. The fact 
that dense cores possess substantial non-thermal internal 
velocities and, in some cases, have values of (3 somewhat be- 
low unity supports the possibility that wave scattering on 
a dense core may lead to substructure like that of TMC-1 's 
Core D. 

In the present paper we examine the effect of MHD 
waves on dense cores. Additionally, we derive the proper- 
ties of the waves that produce substructures within dense 
cores. The paper is structured as follows. In Sect. [2] we de- 
scribe the numerical model. We then describe how the dense 
core evolves in the absence of MHD waves (Sect. [3} in order 
to understand the effect of the waves (Sect.HJ). Finally, we 
summarise and discuss these results in Sect. [5] 



2 NUMERICAL MODEL 

We can use different approaches to study the effect of MHD 
waves on a dense core. One way is to continuously excite 
a wave at the boundary of the computational domain. The 
wave then propagates to the core and interacts with it. An- 
other option is to embed a dense core in a wave. In our 
calculations we use the latter as dense cores arise where the 
waves are. 

An initially uniform quiescent medium is perturbed by 
a magnet osonic wave. We assume that the background mag- 
netic field is in the x — y plane and that wave is propagating 
in the positive x-direction. Furthermore, we work in units 
so that the unperturbed background state is given by 

p = 1, v = 0, B x = I, B y = a. 

The gas pressure p g is given by an isothermal equation of 



state, i.e. p g = pa 2 where a is the isothermal sound speed. 
In our calculations, we use a = 0.1 unless otherwise stated. 

We superpose a fast-mode magnetosonic wave on the 
background state. We do not examine slow-mode waves, as 
they only generate high-density contrasts if their initial den- 
sity perturbation is already large. Contrary to FH02 and 
VFH06, we do not only consider small- amplitude waves but 
also finite-amplitude waves. Then it is no longer possible to 
calculate the initial state using the linear approximation of 
the wave. Therefore, we use the approach of LFH05. For a 
fast- mode wave, the form of the wave written in terms of 
the primitive variables 



P = [p, v x 
must satisfy 
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where 77 is the right eigenvector for a fast wave propagating 
in the positive x-direction, Cf the fast magnetosonic speed 
and A/ = pc 2 — B 2 . Given the x-dependence of one of the 
primitive variables, the others can be determined from the 
above expression. 

As the x-component of the velocity disappears when 
By = 0, it is convenient to specify the profile of v y (x). We 
assume that the 2/- component of the velocity associated with 
the wave varies sinusoidally with x and with wavelength A, 
i.e. 



An, 



2tt(x — x{) 
A 



where A y is a constant and x\ the left boundary of the com- 
putational domain. A y is chosen so that the maximum value 
for the total velocity (\/ v 2 + Vy) is equal to a given ampli- 
tude A. The initial state at a position x + Ax and y can 
then be calculated using 



P(x + Ax,y,0) = P(x,y,0) 



TfAf dv y 

CfB x B y dx 



(x)Ax, 



where the right hand side is evaluated at x. 

We still need to embed a dense core in the fast-mode 
wave. We represent the core by a uniformly dense region 
of radius r at the centre of the computational box. The 
density po of the core is taken such that f3 within the 
core i s between 0.1 and 1 as polar i sation observations sug - 
gest (|Ward -Thompson et al ] l2000l : IWard-Thompsonl l2Q02h 
for some dense cores. As (3 = 2pa 2 /B 2 , the core density 
must lie between 5 and 50 for a = 0.1. Note that there is 
a small variation of f3 within a dense core as the magnetic 
pressure changes across the core. 

Our model can then be specified by the isothermal 
sound speed a, the value a for background B y field, the 
wave properties (i.e. the wavelength A and the amplitude 
A), and the core properties (i.e. the density po and radius 
r). However, as we use dimensionless units, the radius of the 
core and the wavelength of the fast- mode wave are not in- 
dependent parameters. The relevant free parameter is A/r. 
Therefore, we assume a fixed core radius of r = 2 in all our 
calculations and only vary A. 

To solve the 2D MHD equations we use an adaptive 
mesh refinement code. The basic algorithm is a second or- 



Dense core compression and fragmentation 3 



der Godun ov scheme which uses a linear Riemann solver 
(Falle 1991 ). We i nclud ed the divergence cleaning algorithm 
of iDedner et al ] (|2QQ2h to stabilise the numerical scheme. 
The code uses a hierarchy of grids such that the grid spac- 
ing of level n is Ax/2 n where Ax the grid spacing of the 
coarsest level. A refinement criterion determines where in 
the computational domain a finer grid is needed. The com- 
putational domain is — 10 < x < — 10 and — 10 < y < 10 and 
we use periodic boundary conditions. The finest grid has a 
resolution of 800x800. 



3 EVOLUTION OF A DENSE CORE IN THE 
ABSENCE OF MHD WAVES 

In an isothermal gas a dense core is overpressured, which 
means that the core expands. As this introduces waves and 
shocks into the background medium even when the core is 
not embedded in a wave, it is convenient to discuss this ex- 
pansion in the absence of an embedding wave. This problem 
is the magnetic version of the 2D isothermal Riemann prob- 
lem. 

In a purely hydrodynamical case, a shock develops and 
propagates outward into the rarefied gas. At the same time, 
a rarefaction wave reduces the density inside the core by 
travelling inward. As there is no preferred direction, the ex- 
pansion is isotropic. However, this is no longer true in a 
magnetised medium because waves (and, thus, also shocks) 
have a propagation speed which depends on the angle be- 
tween the magnetic field and the direction of propagation. 
Figure [1] shows the expansion of a dense core in a magnetic 
field. 

Slow-mode waves or Alfven waves cannot propagate 
perpendicular to the magnetic field. Thus, only fast-mode 
features develop in those directions. A weak fast-mode shock 
is propagating outward into the background medium, while 
a fast-mode rarefaction wave propagates into the dense core. 
The magnetic pressure increases in material through which 
the shock and rarefaction wave have passed. A tangential 
discontinuity exists between the material that has passed 
through the front of the rarefaction wave and the mate- 
rial that has passed through the shock. The total pressure 
(p g + B 2 /2) is continuous at the tangential discontinuity but 
neither p g or B 2 /2 is. The fast-mode shock is sufficiently 
weak that it does not produce a large density discontinuity. 
Hence, the largest jump in density is at the tangential dis- 
continuity. As it coincides with the edge of the expanding 
dense core, the tangential discontinuity moves outward with 
the gas swept up by the fast-mode shock. 

Slow-mode waves do exist for directions oblique and 
parallel to the magnetic field. A slow-mode shock moves out- 
ward and a slow- mode rarefaction wave propagates inward. 
These slow-mode features lie between the fast-mode shock 
and rarefaction wave. Although the tangential discontinuity 
disappears, there is still a sheet in between the slow-mode 
shock and rarefaction wave separating the initial core mate- 
rial and the shocked background material. 

The density of the material within a region bounded 

1 A discontinuity is a tangential discontinuity if there is neither 
magnetic flux nor mass flux across it, i.e. u n = and B n = 




(Burgess 1995). 



Figure 1. The density structure of an expanding dense core. The 
expansion depends strongly on the angle between the propagation 
direction and the magnetic field. Perpendicular to the field a fast- 
mode shock (FS) propagates outward and a fast-mode rarefaction 
(FR) inward. These waves are separated by a tangential discon- 
tinuity (TD). For expansion at oblique angles the TD breaks up 
into a slow-mode shock (SS) and a slow-mode rarefaction (SR). 
Parallel to the field the FS disappears. The model parameters are 
po = 10, a = 0.316 and a = 0. The core is not embedded in a 
wave. 



by the slow-mode shock is comparable to that of the dense 
core material that has passed through the slow-mode rar- 
efaction wave, which would make the observational identi- 
fication of the dense-core boundary difficult. Since the slow 
magnetosonic speed increases when the angle with the mag- 
netic field decreases, the expansion of the core is faster along 
the magnetic field lines than perpendicular to it. 

In our calculation the fast-mode shock weakens when 
the angle between its direction of propagation and the mag- 
netic field decreases. For parallel expansion the fast-mode 
shock even turns into a linearised wave. Also, the slow-mode 
rarefaction wave merges with the fast-mode rarefaction (see 
Fig. P. 

Finally, the core disperses as there is no force holding it 
together. The dispersal is faster for higher initial densities 
of the core. Also, when the thermal gas pressure becomes 
more dominant in the background gas (i.e. f3 increases), the 
dense core disappears more quickly. 



4 A DENSE CORE EMBEDDED IN AN MHD 
WAVE 

When a dense core is embedded in a fast-mode wave, density 
substructures can be generated by the interaction between 
the wave and the core. The analysis in FH02 shows that 
high- density contrasts are associated with slow-mode waves 
when /3 is small. They also showed that non-linear steepen- 
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Figure 2. The temporal evolution of the maximum density 
within the dense core. The thin solid line shows the evolution 
in a homogeneous background gas (i.e. A = 0), while the thick 
solid line represents a model with a small-amplitude wave. The 
model parameters are: po = 10, a = 0.25, A = 20, A = 0.05 and 
a = 0.1. Initially, (3 = 0.2 in the dense core. 



ing of a fast-mode wave can readily excite these slow-mode 
components. 

To produce substructures, slow-mode waves must be ex- 
cited within the core itself (as their propagation speed is 
low). This means that the time-scale for non- linear steepen- 
ing tns ~ X/2A must be smaller or comparable to the Alfven 
crossing time tA = v/ca- We immediately see that for large 
wavelengths t ns ^> tA even if their amplitude is compara- 
ble to the Alfven speed. Short- wavelength waves, however, 
steepen within the core itself as long as their amplitude is 
not too small. 

The interplay between MHD waves and a dense core is, 
thus, wavelength- dependent. Therefore, it is convenient to 
discuss the results for long and short waves separately. We 
typically take A = 20 and A = 2 in our calculations. 



4.1 Results for large wavelengths 

Although the time-scale for non-linear steepening for large- 
wavelength fast-mode waves is too long for exciting slow- 
mode components within the dense core, the wave strongly 
affects the evolution of the core. The most direct effect of the 
core on a fast-mode wave is the reduction of the propagation 
speed within the core because c/ « B / \fp (for a < c a ). This 
means that a part of the wave - i.e. the part within the core 
- moves slower than the rest of the wavefront. Consequently, 
the velocity of the gas must change at the boundary of the 
core. 

However, more importantly, this induces a change in the 
magnetic field at the boundary of the core, which follows 
directly from the induction equation 



- = Vx(vxB). 

More specifically, B x increases at the top region of the core 
and decreases at the bottom (or inversely depending on the 
sign of the velocity changes). The ^/-component of the field, 
By, shows a similar behaviour, but here left and right of the 
core. 

Thus, the interaction with the wave causes an increase 
of the magnetic pressure. While there is initially no dis- 
continuity in the total pressure perpendicular to the mag- 
netic field (see Sect. [3}, the total pressure outside the core 
can be much higher than inside the core. The magnitude 
of this change depends on the amplitude of the fast-mode 
wave. Small-amplitude waves are too weak to produce any 
significant change and the core then disperses as though no 
fast-mode wave is present in the background gas (see Fig. [2}. 
However, when the amplitude of the fast-mode wave is large, 
the magnetic pressure increases significantly. Now the region 
outside the core is overpressured compared to the core itself. 
A fast-mode shock arises and propagates into the core. This 
shock is weak and, therefore, compresses the gas by a small 
factor as can be seen in Fig. [3Jl. 

Although the expected increase in density from the mo- 
tion of the wave across the core is small, Fig. [3Jl shows that 
the maximum density increases significantly for t > 20. We 
also find that most mass is concentrated within a single re- 
gion. As the time-scale corresponds to the time for a wave to 
steepen into a shock (t ns « 20 for A = 20 and A = 0.5), we 
expect that this high- density structure is generated by the 
interaction between the fast-mode shock and the core ex- 
citing slow-mode waves. However, /3 needs to smaller than 
unity for this process to be effective (see the analysis in 
FH02). Figure [3)3 shows that f3 for the fluid element as- 
sociated with the maximum density obeys this constraint 
reasonably well. 

We can now examine how the model parameters influ- 
ence the wave-core interaction. Of all parameters the wave 
amplitude is perhaps the most important. The wave ampli- 
tude determines how much the magnetic pressure changes 
at the boundary. Small-amplitude waves do not affect the 
evolution of the core, while large- amplitude waves modify 
the entire core structure. In the latter case, a core can be 
compressed to a density considerably higher than the initial 
one. The maximum density is limited by the /3 < 1 con- 
straint. Because of this constraint, it is also more difficult 
to produce high-density contrasts in models where po and a 
are high. Then f3 is already close to unity or higher. 

When the wavelength is larger than the radius of the 
core, the initial position of the core centre with respect to a 
wave node becomes an additional free parameter. However, 
the fast-mode wave still propagates with a lower speed in- 
side the core. The variations described above always occur. 
Hence, the results are essentially independent of the initial 
position. 

4.2 Results for small wavelengths 

A small- wavelength wave steepens much faster into a shock 
than one with a large wavelength. If the amplitude of the 
small-wavelength wave is not too small, the time-scale for 
non-linear steepening is shorter than the Alfven crossing 
time. This means that a shock arises before the wave leaves 
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Figure 3. (a) Similar to Fig. [2] but now for a large-amplitude 
fast wave, (b) The plasma beta associated with the fluid element 
that has the highest density. The model parameters are given by 
po = 10, a = 0.25, A = 20, A = 0.5 and a = 0.1. 



the dense core. Then the generation of density substruc- 
tures can be readily described with the models of FH02 and 
VFH06. 

When a fast- mode wave steepens into a shock, slow- 
mode waves are excited. This produces persistent pertur- 
bations with large density contrasts. The fast-mode wave 
is most efficient in exciting the slow-mode components for 
modest values of a, which is the tangent of the angle between 
its direction of propagation and the magnetic field. Steep- 
ening is slower for very small values of a and the shock does 
not produce slow- mode waves, when a is too large (FH02). 



Figure 4. Similar to fig. [3] but now for A = 2. 

The maximum density due to the non-linear steepening 
is roughly 



Pmax — 



A 2 p 



te « ~T • 



Small-amplitude waves do not change the density in the core 
significantly. The dense core then disperses as in the case of 
small-amplitude large- wavelength waves, i.e. as if the core is 
not embedded in a wave. Large-amplitude waves, however, 
generate high-density regions within the core (see Fig. |4|. 
These density perturbations are aligned with the ?/-axis as 
can be seen in Fig. [5] This is because, except close to the 
boundary, there is little variation in the ^/-direction inside 
the core. 

The density perturbations generated by non-linear 
steepening subsequently decay. This happens on a time-scale 
(see FH02) 

A 
A ' 

Although the dense regions expand and disperse, high- 
density regions can arise for some time due to the inter- 
action of the dense regions with the initial fast-mode wave 
(see VFH06). This interaction excites slow-mode waves and 
produces dense substructures. Since the sheet-like structures 
disperse, the core has a more homogeneous density structure 
with small-scale structures nested inside. This can be seen 
in Fig. [6] Also note that, although large density contrasts 
are generated within the core, the core itself still exists as 
an entity. However, because of the fragmentation of the core 
and associated motions, the kinetic energy inside the core 
contributes significantly to the total energy. The dense core, 
therefore, expands more rapidly than when there is no scat- 
tering wave. 

High-density structures do not only arise inside the core, 
but also at the boundary. In the explanation above, we ne- 
glected the effects occurring at the boundary of the core. In 
a similar way as for the large- wavelength waves, overpres- 
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Figure 5. The density structure at t = 3 inside the dense core. 
The high-density regions are produced due to non-linear steep- 
ening of the fast-mode wave. The initial state is given by A = 2, 
A = 0.5, T e = 0.01, a = 0.25 and p = 10. 




Figure 6. Similar to Fig. \E\ but now for t = 18. Dense sub- 
structures form within the expanding core due to the interaction 
between dense regions and the initial fast mode. 



5 SUMMARY AND DISCUSSIONS 

We have studied the interaction between magnetosonic 
waves and dense cores for which the ratio, /3, of thermal 
gas pressure to magnetic pressure is initially 0.2. We follow 
the evolution of a core embedded in a fast-mode wave. We 
find that large-amplitude waves can change the evolution of 
the core considerably. 

The interaction between the fast-mode wave and a 
dense core depends strongly on the relative size of the core 
to the wavelength. When the core radius is smaller than 
the wavelength, the strongest effects are induced near the 
boundary of the core. Eventually, these effects disrupt the 
global structure of the core, but confines it. Furthermore, 
we find that a core can be compressed significantly. This 
result suggests that a large- wavelength fast-mode wave can 
trigger the collapse of a gravitationally unbound dense core, 
which is of relevan c e to t he turbulent cooling flow picture of 
iMvers fc Lazarianl (|l998h . However, additional calculations 
which include self- gravity are required to confirm this. 

Short- wavelength waves, on the other hand, play an im- 
portant role in the generation of substructure in a core with- 
out breaking up the core. However, contrary to the large- 
wavelength waves, the dense core expands faster than with- 
out the scattering wave. The wave thus does not confine the 
core. Although /3 is close to unity in our simulations, slow- 
mode waves excited by the non-linear steepening of the fast- 
mode wave produce high-density contrasts within the core. 
Such small-scale features have been observed in the cores 
of cold dense clouds such as TMC-1 (IPeng et alJll998h . We 
do not find the large number of microclumps inferred from 
observations of Core D, but this can be readily explained 
by the limitations of our model. The most important limi- 
tation is perhaps that we only studied the effect of a single 
wave. A real velocity field is more accurately described as a 
superposition of an ensemble of waves. Due to wave- wave in- 
ter acti^m^J^g^^ break up into smaller ones. 

lHartquist, Williams fe Vitil (|200lh showed that the sub- 
structures in Core D need to be about 0.1 Myr old to comply 
with the observed molecular variety and abundances. If the 
substructures are being generated due to slow-mode exci- 
tation, the timescales should be of the same order. For a 
core radius of R = 0.1 pc and a — 0.3 km s _1 and a some- 
what larger Alfven speed, the substructures in our simula- 
tions arise on timescales of a few 10 5 yr (a few times R/c a 
which is the relevant time-scale in our simulations). As both 
timescales are of the same order, the microstructure in dense 
cold core can indeed be generated by small- wavelength mag- 
netosonic waves. 
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sured regions arise producing fast-mode shocks which move 
into the core. Because of the compression caused by the fast- 
mode shock, dense substructures, thus, not only form inside 
the core but also at the boundary of the core. 
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